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ABSTRACT 

A new representation for solutions of Maxwell's equations is derived. Instead of being 
expanded in plane waves, the solutions are given as linear superpositions of spheri- 
cal wavelets dynamically adapted to the Maxwell field and well-localized in space at 
the initial time. The wavelet representation of a solution is analogous to its Fourier 
representation, but has the advantage of being local. It is closely related to the rela- 
tivistic coherent-state representations for the Klein-Gordon and Dirac fields developed 
in earlier work. 
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§1. Introduction 

In previous work [1,2] coherent-state representations have been developed for various 
relativistic systems, including Klein-Gordon and Dirac fields and their wave functions. 
The main tool was the analytic-signal transform (AST) [3,4], which gives a canonical 
extension of fields from real spacetime to complex spacetime. When the field f{x) is 
free, its extension f{z) is analytic in the double tube domain 

r = {z = (z, zo) = x-iyeC^\y^ = yl - > 0} = r+ U r_, (1) 

where ±j/o > in 7±. The restriction of f{z) to the forward tube {yo > 0) then 
contains only the positive-frequency part of /, while its restriction to the backward tube 
T_ {yo < 0) contains only the negative-frequency part. The points of T parametrize a 
system of relativistic coherent states e^, with z E and z E T_ representing particles 
and antiparticles, respectively, and it was shown that the state labeled hy z = x ~iy 
has an expected position x at time xq and an expected energy-momentum proportional 
to y. Certain six-dimenional submanifolds of T can therefore be interpreted as phase 
spaces. Furthermore, when the field has a positive mass, it can be reconstructed from 
its values on any one of these phase spaces. However, this reconstruction fails when the 
mass vanishes, due to the divergence of the relevant integrals. (The massless representa- 



tion is not square-integrable with respect to the 'Liouville measure' on phase space.) In 
this paper we develop a square-integrable representation for the most common massless 
system, namely the electromagnetic field. This will be done by reducing the dimenion- 
ality of the "phase space" from six to four. Namely, an electromagnetic field will be 
reconstructed from its values in the Euclidean region (real space- and imaginary time- 
coordinates), as obtained by applying the AST to the field in real spacetime. The result 
will be seen to be a multidimensional generalization of wavelet analysis. 

That relativistic coherent states behave like wavelets due to Lorentz contractions 
has been noted in [3]. However, only in the massless case can the correspondence be 
complete, since a positive mass provides a scale, namely the width of in its rest 
frame. This led us to the expectation that massless fields (where no rest frames and 
canonical scales occur) possess a natural wavelet representation [4]. The simplest case 
is the wave equation in two spacetime dimensions, for which a wavelet representation 
was indeed constructed in [5]. As shown there, the wavelets are covariant not only 
under the group of affine transformations (translations and dilations) but also under 
the larger group of conformal transformations. Here we develop the simplest aspects 
of a similar construction for the electromagnetic field in four spacetime dimensions. A 
more detailed analysis of this and related topics, such as conformal invariance and the 
construction of polarized electromagnetic wavelets, will be dealt with in a forthcoming 
paper [6]. To the author's knowledge, the results given in [5] and in the present paper 
represent the first successful "relativistic coherent-state" formulations of massless fields 
(in two and more spacetime dimensions, respectively). 

It is remarkable that a single function, namely the AST-extended field /{z)., com- 
bines the concept of phase space (parametrized by position and momentum) with the 
concept of wavelet space (parametrized by position and scale). This unification is ac- 
tually the "relativistic dividend" earned from the construction of relativistic coherent 
states, in the same way as relativity unifies energy with momentum and electric fields 
with magnetic fields. To see this, note that although y is proportional to the expected 
energy-momentum of e^, it cannot be an energy-momentum since it has units of length 
(which scale oppositely to those of momentum). Rather, y/yo gives the expected veloc- 
ity, while the imaginary time yo gives a scale parameter. In the Euclidean region, y = 0, 
which means that represents a spherical wave which first implodes towards a point 
in space, then explodes away from it. Then yo measures the scale of this wave by giving 
its diameter at the instant of maximal localization. Such e^'s will form our system of 
electromagnetic wavelets in the next section, and the field itself will be constructed from 
them. The relation of our formalism to the usual one-dimensional wavelet analysis is 
discussed in Section 3. 

We consider solutions of Maxwell's equations in free spacetime R^. It will be 
convenient to unify the electric and magnetic fields E(a;) and B(a;) into a single complex 
vector field F(x, t) = F{x) = E(a;) +i'B{x), which then satisfies (with the speed of light 
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c=l) 



V • F = 0, idtF = V X F. (2) 

Let us review the usual solution using Fourier transforms. The above equations imply- 
that F satisfies the wave equation dt ^F = V^F, hence it has the form 

F{x) = f dpe-'P^'fip), (3) 
Jc 

where = Po — P^, px = pot — p • x, C = C+ U C_ is the double light cone with 
C± = {{p,Po) I ±Po = |p| > 0} (the origin p = is excluded) and dp = {27r)~^d^p /2|p| 
is the Lorentz-invariant measure on C. Maxwell's equations imply p ■ f{p) — and 
Pof(p) = ip X f(p) for p in C, which can be solved by introducing a four-potential. 
A Poincare-invariant norm (and associated inner product) on solutions is given [7] in 
momentum space by 

||ff ^ / dpu;-^\f{pr, (4) 
Jc 

where to = |p|, and we denote the Hilbert space of all solutions with finite norm by 
Ti. Gross [7] has shown that a norm unitarily equivalent to the above is also invariant 
under the fifteen-dimensional conformal group C of spacetime, which is generated by 
the Poincare group together with uniform dilations (x — > ax, a > 0) and inversions 
in the unit hyperboloid {x — > x/x^). This gives a unitary representation of C on H. 
(That Maxwell's equations are invariant under C has been known for a long time [8]; 
the unitary representation gives that invariance a quantum-mechanical flavor, since the 
elements of 7i may now be regarded as single-photon wave functions.) The new norm 
has the following non-local expression in terms of the fields at time t = 0: 

^ Jnp \^ y\ 

where the asterisk denotes complex conjugation. In the next section we shall find a 
simpler expression for ||F|p in terms of the values F(x, —is) of the field in Euclidean 
spacetime (real space but imaginary time). This will lead us directly to the wavelet 
expansion of F. 

§2. The Wavelet Representation of Solutions 

Our extension of F(x, t) to complex spacetime employs the analytic-signal transform, 
developed in [1,2,3] and further applied in [4,5]. Given an arbitrary but reasonable 
function /(x) on R"^ (a smooth function with mild decay is more than sufficiently 
reasonable), its analytic-signal transform is the function /(x + iy) on C"' defined by 
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1 poo J 

/(x + zy) = -/ ^/(x + ry). (6) 
In terms of the Fourier transform / of /, 

/(x + zy) = {2ti)-^ j d-p2^(p-y)e^P-(-+^y)/(p), (7) 

where ^ is the step function {d{u) = 1 if m > 0, ^(0) = 1/2 and 9{u) = if tt < 0). 
Although /(x + iy) is in general not analytic, it can be easily shown to be partially 
analytic in the direction of y 7^ 0. Furthermore, if /(p) has certain support properties 
(vanishes outside of a solid double cone V C R""), then /(x + iy) is analytic in a 
corresponding domain in C"^ (the double tube over the double cone V dual to V). This 
will be seen explicitly below, where R" is spacetime. 

Application of the AST to the electromagnetic field F(a;) gives 

F{x -iy) = — / F{x -Ty)= / dp2e{j)y) e^^f^^-^^^) f(p). (8) 

TTZ T-l Jc 

(We write x — iy instead of x + iy to conform with the convention used in Streater and 
Wightman [9].) If y belongs to the forward solid light cone VI (i.e.,y2 >Oandyo >0), 
then 9{py) = on C_ and 0{py) = 1 on C-|-, hence F{x — iy) is analytic at x — iy 
and contains only the positive-frequency part of the field. Similarly, if y belongs to the 
backward solid light cone V!_, then F{x — iy) is analytic at x — iy and contains only the 
negative-frequency part of the field. Thus F{x — iy) is analytic in the double tube T 
defined in Eq. (1). The jump discontinuity in F(a; — iy) between the future and past 
imaginary directions is related to the multidimensional Hilbert transform of F{x) (see 
[4,5] for details). In this paper, we shall be interested only in the transform of F(x, t) 
to complex time t ^ t — is (i.e., F{x — iy) with y = (0, s)): 



F(x,t-is)= / dp2e{pos) e-'P°^^-''^+'P-'' i{p) 
'c 



Jr^ (27r 



(27r)^a; 



e{s) e-"(^+^*) f (p, iu) + 0{-s) e"(^+^*) f (p, -to)] (9) 



= |a;-^ [e{s) e-'^(^+^*) f (p, u) + e{-s) e'^(^+**) f (p, -w)] } --(x), 

where " denotes the inverse Fourier transform in R^. Note that F(x, t — is) is analytic 
in t — is whenever s 7^ 0. The last line in Eq. (9) implies, by Plancherel's theorem, that 
for s^O, 
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/.3x|F(x,*-i.)P 



Irs {2nfu;^ 

= / ^ 

Jus (27r)^cc;2 



e{s) e-'^(^+^*) f (p,a;) + e{-s) e^^'+'^'^ f (p, -a;) 
{e{s) |f (p, + e{-s) e^- |f (p, -a;)^} , 



(10) 



and therefore 



/ d^:x.ds\F{x,t - is) 



R3 (27r)"'2c^3 

/ dp uj-^ \i{p)\^ 
Jc 



{|f(p,a;)|2 + |f(p,-a;)|2} 
llf||2. 



(11) 



Hence we define 



IFP = / cZ^xcZs|F(x,t-is) 

'R4 



(12) 



so that ||F|| = ||f||. Note that the new norm ||F|| and its associated inner product 
(G|F) are local in the Euchdean spacetime variables (x, s). Define the function 
^y,t-is{p) on C by its complex-conjugate as 



The spacetime function corresponding to By -is is 



ey_is(x,t)= dpe *^^ey_is(p) 
Jc 



(13) 



/ #a;2 2d{pos) e-Po(3+^^)+^P•(x-y)_ 



(14) 



This is a scalar solution of the wave equation which depends only on |x — y | and t — is, 
being analytic in the latter variable whenever s ^ 0. Although Cy^t-is does not belong 
to H (since it has no polarization, being a scalar), we shall write Eq. (9) as 



F(x, t-is) = { e^,t-is I f ) = ( ex,t-is | F ), (15) 

where ( | ) denotes the inner product in Ti, expressed either in momentum space or, 
equivalently, in Euclidean spacetime. Then it follows from ||F|p = ||f|p that the inner 
product of two solutions G and F in 7i can be written as 
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(G|F)= / d^ydsG{y,is)*-F{y,is) 

7 (16) 

which gives the continuous resolution of unity 

/ d^yds \ey-is){ey-is\ = In, (17) 

where is the identity operator in the space of solutions and the equahty holds in a 
weak sense. Hence for F in we have 

F(x, t - za-) = ( ex,t-ia I F ) = / d^yds{e^,t-ia\ey-is){ey-ia\F). (18) 
A straightforward computation gives 

{e^,t-ia\ey,-is) = ^2 (^2^^2)3 ' T = S + (7 + zt, r = |x - y | . (19) 

To obtain the wavelet expansion of F(x, t), we simply take cr = in Eq. (18): 

F(x,t)= / d^yds{e^,t\ey-is)F(y,-is) 

7 (20) 
= / d^ydsey-is{x,t)F{y,-is). 

As we have seen, Cy _is(x, t) satisfies the wave equation in (x, t). Since it depends only 
on r = |x — y|, it is a spherical wave centered at x = y. Furthermore, s acts as a scale 
parameter, since 

ey -is(x, t) = Cy/s _j(x/s, t/s). (21) 

Hence it suffices to examine any one of the wavelets Cy^-is- Figures 1-4 show the behavior 
of the wavelet centered at the origin with scale s = — 1, i.e. of the basic wavelet 

w{r,t)^n'eo,i{^,t) = ^^0^, r^|x|. (22) 

These figures confirm that for t < 0, ey,_is(x, t) is an incoming spherical wave which 
builds up rapidly (at the speed of light!) to a well-localized packet in the ball |x — y | < 
\/3 \ s\ and decays rapidly for t > into an outgoing spherical wave. The characteristics 
t = ±r appear as ripples in figures 1-3. These properties partly justify our use of the 
term "wavelets." Further justification is given in the next section. 
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3. Relation to Wavelet Analysis 

The correspondence F(x, t) F(y, is) is analogous to the Fourier transform F(x, t) — » 
f(p), and the reconstruction of F(x, t) from F{y, is) is analogous to its reconstruction 
from i{p) via the inverse Fourier transform. The role of the plane waves e*^^ is now 
played by ey^_js, and the local nature of these functions means that the behavior of 
F(y, zs) is correlated with that of the field F(x, t) in real spacetime, the correlation 
being strongest when t = since Cy^-jg is then most localized. [By contrast, the Fourier 
transform has no such local property since the plane waves extend to all of spacetime; 
thus a small local perturbation in F(x, 0) can cause an unrecognizable change in f (p).] 
In fact, F(y, zs) is a version of F(x, 0), blurred to resolution \s\. Eq. (20) with t = 
states that F(x, 0) is recovered by superposing wavelets with different centers y and 
scales s, with F(y, is) as the coefficient function. This representation of functions as 
superpositions of their blurred versions is typical of wavelet analysis [10,11,12] and the 
related multiresolution analysis [13,14]. Moreover, in our case the wavelet representation 
also gives the time evolution because our wavelets are dedicated [4,5] to the dynamics 
of the electromagnetic field. 

In momentum space, the wavelet with s — 1 and y = is represented by 

eo,-iip,po) = 2uH{po) e-^'o, p e C. (23) 

This is a multi-dimensional generalization of a basic wavelet which has appeared pre- 
viously in the literature in connection with the usual (one-dimensional) wavelet theory 
[15,16], and also in a wavelet analysis of solutions of the wave equation in two space- 
time dimensions (Ref. [5], Eq. 161). In fact, the AST is a special case of a windowed 
Radon transform, which in turn has been shown to be a multivariate generalization of 
the wavelet transform. 
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